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We study the fermion pair production from a strong electric field in boost-invariant coordinates 
in (3+1) dimensions and exploit the cylindrical symmetry of the problem. This problem has been 
used previously as a toy model for populating the central-rapidity region of a heavy-ion collision 
(when we can replace the electric by a chromoelectric field). We derive and solve the renormalized 
equations for the dynamics of the mean electric field and current of the produced particles, when the 
field is taken to be a function only of the fluid proper time r = y/t"^ — z^. We determine the proper- 
time evolution of the comoving energy density and pressure of the ensuing plasma and the time 
evolution of suitable interpolating number operators. We find that unlike in (l-fl) dimensions, the 
energy density e closely follows the longitudinal pressure. The transverse momentum distribution of 
fermion pairs at large momentum is quite different and larger than that expected from the constant 
field result. 
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I. INTRODUCTION 



The "Schwinger mechanism" for pair production has 
been used in various phenomenological models for par- 
ticle production following a high-energy heavy-ion colli- 
sion. One theoretical picture of high-energy heavy-ion 
collisions begins with the creation of a flux tube contain- 
ing a strong color electric field [T]. The field energy is 
converted into particles such as qq pairs and gluons by 
the Schwinger mechanism [21 [31 H] . This mechanism has 
been implemented in a phenomenological fashion in event 
generators for particle production such as the Lund string 
model of hadronization [5] or the Hijing model [6]. More 
recently another picture of heavy-ion collisions, based on 
the color glass condensate model of high density for quan- 
tum chromodynamics (QCD) [71 [S] has been put forward. 
This model leads to the picture that a heavy-ion collision 
produces an initial semi-classical chromoelectric field in 
the longitudinal direction. Kharzeev tt al. [S] have shown 
that if one looks at a perturbative parton cascade model 
and studies inclusive production of gluons in a gluon cas- 
cade, that this is equivalent to the production of a gluon 
from a background classical chromoelectric field in the 
longitudinal direction. This recent work gives credence 
to the idea that as far as gluon production is concerned, 
one can replace the dynamics of heavy-ion collisions by an 
initial condition on a semiclassical chromoelectric field. 
In these recent papers however, no attempt has been 
made to actually study the time evolution of the result- 
ing plasma and the backreaction of the production on 



the initial chromoelectric field. Some early studies had 
been done phcnomenologically on this type of problem 
using a kinetic theory model in which a relativistic Boltz- 
maim equation is coupled to a simple Schwinger source 
term |10l [TTl [T^ [T?j . and a Wigner function transport 
approach for an SU(2) version of QCD was recently done 
by Skokov and Levai |14j . A first principle (quantum field 
theoretical) calculation for pair production and backreac- 
tion from strong fields was done by one of us and collab- 
orators in the appropriate kinematics for heavy-ion col- 
lisions in (1+1) dimensions in an abelian approximation 
where one ignored the color degrees of freedom. The rea- 
son for revisiting this problem now is two-fold. First, an- 
alytic results for the transverse momentum distribution 
functions for particles produced by constant electric and 
chromoelectric fields have recently been obtained [TSlllSj . 
For the constant chromoelectric field, the results for pair 
production are different than for an electric field in that 
the transverse distribution of jets depends not only on 
the energy density of the field but also on the direction 
the field is pointing in color space, i.e. the color hy- 
percharge. Thus it is important to know first, how the 
backreaction affects the transverse momentum distribu- 
tion function both for quantum electrodynamics (QED) 
and QCD and secondly whether adding interactions in 
a 2-PI 1/N expansion will modify the one loop result. 
Here we will address the problem of finding the trans- 
verse distribution function for the abelian case in (3+1) 
dimensions in a realistic kinematic scenario. The QCD 
problem will be addressed in a separate paper. 
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First let us review the history of both analytic ap- 
proaches to the constant electric field and chromoelectric 
field problem as well as numerical studies of the back- 
reaction problem. In his 1951 classic paper, Schwinger 
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derived the following one-loop nonperturbative formula 



dW 



47^3 



E 



-mrm^ /\eE\ 



(1.1) 



for the probability of fermion pair production per unit 
time per unit volume from a constant electric field E via 
vacuum polarization f?^ by using a proper time method. 
The result of Schwinger was extended to QCD by Claud- 
son, Yildiz and Cox [17]. However the pt distribution of 
the 6+ (or e~) production, dW/d'^xd'^pr, could not be 
obtained using the proper time method of Schwinger. A 
WKB approximate method was used for this purpose by 
Casher et al. ^18, , but an exact method to do this problem 
(of determining the transverse distribution of pairs) was 
not found until recently [T5l [B]. For QED the WKB 
analysis gave the correct answer which depended only 
on the energy density of the electric field. However, for 
QCD, the WKB answer was incorrect for QCD in that it 
did not contain the second Casimir invariant of 5J7(3), 
C2 = [dabcE''E^E''f, as shown in Refs. HH and US In 
the case of fermions in QED one finds for the transverse 
distribution of fermion pairs: 
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(1.2) 



The purpose of this paper is to consider the backreac- 
tion problem in (3+1) dimensions in a situation which is 
related to the kinematics of particle production by strong 
chromoelectric fields, namely initial conditions where the 
center-of-mass energy is so high that all distribution func- 
tions are boost invariant in the longitudinal direction so 
that physical quantities only depend on the longitudinal 
proper time which is the same in the boost-invariant limit 
to the (l-l-l)-dimensional fiuid proper time. Our goal is 
to see how the original result of Casher et al. [TB] , which 
has been recently rigorously derived by Nayak for the 
transverse distribution of fermion pairs [19 , is modified 
by the expansion of the ensuing plasma and the backre- 
action on the electric field. 

The backreaction problem was first studied numeri- 
cally in real time for both scalar QED and QED by 
Cooper, Mottola and collaborators in (1-1-1) dimen- 
sions [2ni mi [22 and then also in boost invariant co- 
ordinates relevant to heavy-ion collisions in Ref. [15] • 
Ref. |23J, a strong abelian field was used as a model for 
particle production in the central-rapidity distribution. 
In that work, the boost invariance of the problem was 
used to show that many features of the hydrodynami- 
cal model were appearing even though there were no in- 
teractions kept that would lead to equilibration. Also, 
in that paper, although the theory was formulated in 
(3+1) dimensions, numerical results were only presented 
for (1+1) dimensions, so that the transverse distribution 
of secondaries was not studied. To remedy this particular 
deficiency of our previous work, here we investigate the 
dynamics of the particle production as a function of time 



in cylindrically-symmetric boost-invariant coordinates in 
(3+l)-dimensional QED. 

The present paper builds on our previous papers on 
fermion pair production in (1+1) dimensions by strong 
electric fields with a backreaction of the current on the 
field |23l 124] . As in our previous work, we employ quan- 
tum field theory methods in the large-N approximation to 
find the particle production rate. The next logical steps 
are to extend this result to QCD in this one-loop approx- 
imation and then to do a self-consistent resummed 1/N 
expansion to study the competition between the thermal- 
ization of the plasma and the expansion. In this way we 
will build up gradually the machinery to ask important 
questions about the thermalization and expansion of the 
quark-gluon plasma in a model based on the Schwinger 
mechanism. 

The paper is organized as follows: In Sec. [IT] we derive 
the equations needed for this calculation. In Sec. Ill we 



derive the components of the energy-momentum tensor 
in this coordinate system and show that it is conserved. 
In Sec. |IV| we introduce the concept of the quasiparticle 
phase space distribution function. This quantity can be 
extracted from the field theory energy density and then 
used to determine the distributions of pairs produced in 
the center-of-mass frame. In Sec. [V] we discuss our nu- 
merical approach and present results of our calculations. 
We conclude in Sec. |VI[ In App. [A] we explain the nota- 
tion we use throughout this paper. In App. [B] we derive 
the transverse helicity eigenvectors we use in the main 
text of the paper to expand the fermi field, whereas in 
App. [C] we derive an adiabatic expansion of the Dirac 
equation which are used throughout the paper to study 
the large momentum behavior of integrands. 



II. THEORY AND NOTATION 

In cartesian coordinates = (t,x,y,z), and using the 
matric r]ab = diag( 1,-1,— 1,-1), the lagrangian density 
for this problem is given by 



C = m {I'^i^da-e Aa (0 ] - M } ^-(0 

_ 1 

4 



^ F'^\OFab{^). (2.1) 



where F''''(C) = A'' {C) - A" {C) ■ Equations of motion 
are given by 

{ 7"^ [ Ida ~ 9 Aa{i) ]-M) ^(0 = , (2.2) 
daF'^\S,) = {f{i) ) = 5 ( [ m.l' m ] >/2 . (2.3) 

In this semiclassical approximation, we quantize the 
Dirac field while the electromagnetic field is treated clas- 
sically. This approximation can be made precise by con- 
sidering N flavors of quarks interacting with the electro- 
magnetic field and considering the limit where N 00 
after appropriate re-scalings. Systematic corrections are 
given by the 1/A^ expansion as discussed in Ref. i2^ and 
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references therein. We use the standard representation 
of the 7-matrices 



1 
-II ' 



a' 

-a' 



(2.4) 



where cr* are the usual Pauh matrices, and 1 is the unit 
(2 X 2) matrix. 



A. Dirac's equation in boost-invariant coordinates 

Boost-invariant coordinates = {t, p,9 ,rj) are de- 
fined by 

i = TCOsh77, z^rsinhyy, (2-5) 
X — p cos , y = P sin 9 . 

We use Roman indices to indicate the cartesian frame 
and Greek indices for the cyhndrical-hyperboHc frame. 
The connection between the cartesian frame (d^°), and 
the boost-invariant frame (da;^) is described by a vierbein 
matrix which for our case is 

dC = V%{x) dx^ , = V\{x) da , (2.6) 

with 

(cosh 77 T sinh 77 

cos 61 -psin6i 

sin 6* pcos6l 
sinh 77 r cosh 77/ 

The inverse vierbein matrix, which we write as V^a{x), 
is given by 



da;^ = l/^a(x)dr, 5a = y^(x)9^ 



(2.7) 



with 



dxt" 
cosh rj 



(2.8) 



— sinh 77 ^ 
cos 6* sin 6* 
—smO/p cos9/p 

V— sinh77/T cosh7;/Ty 

The metric g^v{x) in boost-invariant coordinates is 



gt.Ax)=ilabV\{x)V\{x) (2.9) 
= diag(l,-l, V,-^')- 

We raise and lower Latin indices by the 77-metric and 
Greek indices by the g-metric. So Dirac's equation in 
boost-invariant coordinates can be written as 

{ 7'^(x) [i 9^ -eA,Xx) ] - M } 7/;(x) = , (2.10) 



where we have defined 7'' (a;) = "i°'V^a{x). In this coor- 
dinate system, the 7-matrices are given by 



7^(2;) = cosh 7/ 7*^ — sinh 77 7'^ , 
7''(a;) = cos 9^"^ + sin 9 7^ , 
7*^(2;) = (-sin6'7^ + cos6'7^)/p, 
7'' (a;) — ( — sinh t; 7° -|- cosh r/j^) /i 



(2.11) 



The fermi field ^^{x) in boost-invariant coordinates obeys 
the anticommutation relation 



{i;,{T,p,9,r^),tPl{T,p',9',r^')} 
VPP 



9')fc^. (2.12) 



It is simpler, however, to solve Dirac's equation in a 
Lorentz-transformed frame which diagonalizes the vier- 
bein. The Lorentz transformation that does this is 



A\{9,7j)^V\{T,p,9,r^)V'^,{T,p) 

/cosh?7 sinh 77 N 

cos 6* -sin 6* 

sin 6* cos 6* 

\sinh7; cosh 777 



(2.13) 



where V^b^r, p) are the diagonal vierbeins given by 

V%{t, p) = V\{t, p, 0, 0) = diag( 1, 1, p, r ) (2.14) 
V\{t, p) EE y^(T, p, 0, 0) = diag( 1, 1, 1/p, 1/r ) . 

We define 7-matrices in this frame with a bar 

f'(r,p) = t/^(T,p) 7'^ = 7''(t,p,0,0). 

They are given explicitly by: 

7- = 7O , 7''(r) = 7Vr , 

f{p)^lVp- 



(2.15) 
(2.16) 



Now let S{9, rj) be an operator which induces this Lorentz 
transformation on the 7"^ matrices in the orthogonal 
frame. 



S-\9,r^)rS{9,r^)=A\{9,r^)j' 



(2.17) 



Then it is easy to show that S{9, 77) = Sp{9) Sp{ri) is given 
by a product of operators, where 



Sp{9)^eM0tYm 

= cos(6'/2) + 7I72 sin(6l/2) , 

^,(77)=exp[77 7V/2] 

= cosh(77/2) + 7°7^ sinh(77/2) 



(2.18a) 
(2.18b) 



Furthermore, from Eq. ( 2.13[ ), we see that S{9,7j) trans- 
forms the ^'^{t, p,9,ri) matrices into the 7^(r, p). 



S-\9,v)r{x)S{9,r^)=r{r,p) 



(2.19) 
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Now let us note that 



B. Mode expansion 



s-\e,ij)~-i^{x)d^s{e,i) 

= 7'^(T,p)[a^ + n^(0,77)] , (2.20) 

where we have defined a connection 11^(0, 77) by 



(2.21) 



The only nonzero connections are when fi = 9 and n — rj. 
So using Eq. (2.181, we find 



ne = 7'7V2, and H,, - 7V/2 , (2.22) 

which are independent of 9 or 77. The covariant derivative 
is given by 

= + n^(x) = d^+ U^ix) + le A^{x) . (2.23) 

ChristofFel symbols for boost-invariant coordinates, 
which we will need later, are given by 



(2.24) 



from which we find the only non-vanishing elements to 
be 



'^\r] ^ I ^Ti-j — -TJjr — J (2.25) 

pB — ^ Bp 



So Dirac's equation (2.10) can be transformed to the 



boost-invariant frame by defining 

i,{x) = S{9,^) ^{x)/V^., 



(2.26) 



and multiplying the equation through by S ^{9,r/), 
which gives the equation 



[z7^(t,p)V^-M] 



7-0. 



(2.27) 



For our case, we assume that the vector potential is in the 
ry-direction and depends only on r, so we choose Afj_(x) = 
( 0, 0, 0, —A{t) ), which defines A{t) as the negative of the 
covariant component. Then ( |2.27 1 simplifies to 



{ij°dr + ij^[dp + l/(2p) ]+t^^de/p 

+ -f^ [idn + eA{T)]/T~-M}4>{x)^0, (2.28) 

which is the equation we want to solve. Here (p^x) field 
obeys the simpler anticommutation relation 

{ 0„(t,x), <^^,(r, x') } = S^,^, , (2.29) 

where x = {p,9,ri). Our notation is explained in App.[A| 



An expansion of the (f>{x) field in terms of transverse 
helicity eigenstates can be carried out using the separa- 
tion of variables methods explained in Ref.|5S]and further 
discussed App. [B] The expansion is given by 



(x) = ^4'V(''(:.), 



(2.30) 



k,A 



where k = {krpk±,m,h) with 



E 



+00 



dkri 



+00 



fcj^ dfc 
27r 



+ 00 



-E E 

m— — 00 h—-^\ 



(2.31) 



(See App. [A|for our notation.) Here A = ±1 labels initial 
positive and negative energy states, ft, = ±1 labels the 
transverse helicity of the state, and m the value of the 
z-component of the angular momentum operator. The 
time-dependent spinor mode functions 'i)'^^ (x) are given 

by 



(2.32) 



where k = (kri,k±,h) and km — {kj_,m), and where the 
transverse helicity eigenvectors Xfe„,/i(p, 0) are given by 



J(m+l/2) e 

V2 



Jm{k±p) 



(2.33) 



In App. |B]in Eqs. ( |B23| ) and ( |B24| ), we show that they 
are orthogonal and complete. 

The 4'''^\{t) mode functions form a two-dimensional 
spinor, 



which satisfies the equations of motion 
where the Hermitian matrix Hkir) satisfies: 



+M '^ri{T) — ihk± 
TTrjir) + ihk± — M 

kfe(r) • a- , 



(2.34) 



(2.35) 



(2.36) 



with TTq^r) = [ fc,j — eA(r) ]/t the fcmeizc momentum, and 
where 



kfe (r) = TT^ (r) ei + hk_i_ 62 + M 



(2.37) 



We define a density matrix pI^\t) and "polarization" 
vector p''^\t) by 



Pk^ (r) - (r) (r) = [ 1 Pr ^ (r) • <T ] , (2.38) 
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so that from Eq. (2.351, the polarization vector satisfies: 

(2.39) 



We find an adiabatic expansion to second order of the 
polarization vector in App. [C| Since Hkir) in (2.361 is 



Hermitian, the length of the spinors (I)1^\t) is conserved 

5r[0i'^^(r)0i''Hr)] =0. (2.40) 

So if we choose the two spinors labeled by A to be orthog- 
onal at r = tq, then they remain orthogonal for all t. In 
Sec. |II C| below we do this, so we can assume that these 
spinors are orthogonal and complete for all values of r 



(2.41a) 
(2.41b) 



Near t = tq = 1/M where we take A{to) = 0, the 
Hamiltonian (2.361 is approximately independent of r. 



Hkir) « Ho-k, where 
+ 1 



Ho;k = M 



where 



^ ^"-f-UMko...., (2.47) 



ko-fe = fc,, ei + /ifcj_ 62 + 63 



(2.48) 



with fcj^ = k^jM. We write the eigenvalue equation for 
the Hamiltonian Ifo-fc as 



tlQ-.k 0o;fe - ^0;fc <P0;fc 



(2.49) 



with J-^j, = XMujQ-k , ^0:fc = + fc]^ + 1 
and where 



Probability conservation also requires that the polariza- 
tion vector f'"^\t) for both of these solutions to remain 
on the unit sphere for all time t. 



Using the orthogonal relations (B23l and (2.41a I, we 
can invert expansion ( |2.30| to obtain for any time t. 



A 



E 



where our notation is explained in App. [A| Using (2.291 
the mode operators obey the anticommutation relation 



k ' 



(2.43) 



It is traditional to define separate positive and negative 
energy mode operators by setting 



- "k ' 



and 



(2.44) 



We choose our initial state to be the vacuum with no 
particles or anti-particles present. Then 



This means that 



and 



(A)t J(A')i\ _ 



'.^L'M) = -A5,.V"k.k'> 



. 



(2.45) 



(2.46) 



a result we will use in Sec. Ill Dl below. 



C. Initial conditions 

There have been several methods used to set initial 
conditions for the fermion field. We investigated two of 
these methods in Ref.[24land came to the conclusion that 
both methods produce essentially the same results, so we 
choose the simpler "one-field" method here. 



i LOQ-k + l f 1 
2wo;fe VCfe 



(2.50a) 
(2.50b) 



(2.42) with Cfe — {krj + ihk±) /{ojQ-k + !)• We use these eigenval- 



ues for initial values of the spinors (fi^^^ (r) at r = tq 



(2.51) 



which defines what we call positive and negative energy 
solutions of the full Dirac equation. Since the initial 
spinors are orthogonal and complete, the full solutions 
of the Dirac are also orthogonal and complete. The den- 
sity matrix p^^l at tq is given by 



fOik — rO;fc VOik 



where the initial polarization vector Pg'^^ is given by 

= Ako;fc/c^O;fc- (2.53) 

D. Maxwell's equation 

In boost-invariant coordinates. Maxwell's equation 
reads: 



(2.52) 



(2.54) 



where — pr. Now A^(x) — (0,0, 0, —A{t) ), so the 
only non- vanishing elements of the field tensor are: 

Fr^^{x) = ~F^^r{x) = -drA{T) = T E{t) (2.55) 
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This last equation defines what we call the electric field 
E{t) = —drA{T) It. Then using the metric g^^{x) — 
diag( 1,-1, -1/p, -l/r^ ), we get: 

F^.''(t) = -P-^'-'^ir) = -E(t)It , (2.56) 

and Maxwell's equation becomes: 

drE{T) = -J{t) . (2.57) 

Here we have defined a "reduced" current J(t) by: 



e r 



^(r)-^([V(2;),7''(r)V'(a:)]) 



([0t(a;),7V^(x) 



2t 



(2.58) 



Using the field expansion (2.301 and the expectation 



value (2.46) of the mode operators, we find for the re- 



duced current: 



2r 



k,k' A, A' 



= -^EM'^l'^'(^)7V'/'L'W^)] (2.59) 



k.A 



Now since 



(73 
(73 



(2.60) 



and using the fact that (73 Xk h(P^ ^) ^ ^fc -h(P^ 
the relation 



+00 



E xL./.(/''^)Xfe„,,(p,^) 

'I— —00 

+00 

= E + =1, (2.61) 



we find from (2.32 1 that the reduced current can be writ- 
ten as 



•^(^) = -27E^Tr[pi^)(r)(7i] (2.62) 

fc,A 



Here we have used the completeness statement (2.41b) 



to write the current in terms of positive energy solu- 
tions only. In the last line, we changed integration vari- 
ables from fc^ to TT^, using dvr^ = dfc,,/T, and defined 
Pfe^.;i(7r,p r) = Pj.(r). Maxwell's equation (2.57) be- 
comes: 



a.£;(r) = e5]pg^_,(7r„r) 



(2.63) 



Recall that P^^^ ^(tt^jt) is the first component of the 



positive energy polarization vector. Eq. (2.39) with ini 



tial condition (2.53), and Eq. (2.63) need to be solved 



simultaneously for the system dynamics. 

As it stands, the integral for the current in Eq. (2.63) 
diverges. We renormalize it using the adiabatic ex- 
pansion of solutions of the Dirac equation we found in 
App [C| Setting e 
Eq. ([2j63| gives 



1 and substituting (C21al into 



E{r) 



(2.64) 



p 



^-(fci + M2)f 



1 Tiv, 
4 uj^ 



5 TT,, TTjj 

8 w7 



where here iv = [tt^ + fc^ + M^]^/^. The dot refers to a 
derivative with respect to r. So the first term vanishes 
by symmetric integration over tt^. For the other terms, 
we note that 



7r^(r) = [fc^ - eA(r)]/T , 

7i'r;(T) = -7r,,(T)/T + eE{T) , 

nr,{T) = 27r^(T)/r2 - eE{T)/T + eE{T) 



So the only terms which survive in (2.64) are 





eE{T) 



27r 



27r 



(2.65a) 
(2.65b) 
(2.65c) 



(2.66) 



( 4tj' 



57r| 
4cj7 



ei;(T) 
4tj5 



Here we have introduced a cutoff A in the fcj^ integral. 
Carrying out the integrals in (2.66) and moving terms 



proportional to E{t) to the left-hand-side, we find the 
adiabatic expansion of Maxwell's equation to be 



[l + e'{5e')]E{T)^eR[eA{T)] 
where 5e^ is given by 



^ln[A/M], 



(2.67) 



(2.68) 



and i?[eyl(r)] is a finite functional of the product eA(r), 
or derivatives of this quantity. We define the renormal- 
ized charge Cr by 

„2 

(2.69) 



l + e2((5e2) 



Then since eyl(r) = erAr{T), the adiabatic expansion of 
Maxwell's equation (2.67) reduces to 



Er{T) — Cr R[er Ar{T)] , 



(2.70) 



which is now finite. We conclude that we can regularize 
Maxwell's equation by subtracting from the integrand 
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the adiabatic expansion of -P^.^^ hi''^V' ''') ^^'^ addition Here = Afj_{x) and the subscript notation (/x, i 



renormalizing the charge. This gives the equation 

A 



drEir) = 



l + e2(Je2) 



E 



(2.71) 



III. ENERGY-MOMENTUM TENSOR 



means to symmetrize the term. From our results for the 
field tensor in Eq. (2.55) in Sec. II D the field part of the 



energy-momentum tensor is given by 



T^f^ = ^ diag( E\ E\ p'E\ ~t'E' ) 



1 



(3.3) 



We denote the matter part of the energy-momentum ten- 
sor as: 



In the boost-invariant coordinate system, the average 
value of the total energy-momentum tensor is given by 
Eqs. (4.1) and (4.2) of Ref. [231 and is the sum of two 
terms (notice sign convention) 



rp rrimattcr i T-ificld 



where the matter and field contributions are given by 
1 



(3.1) 



-([ ^(x), 7(^(x) {z D,^ ^{x)) ] + h.c. ) (3.2a) 



1 



rpficld _ _ p.ap p J, p a/3 p 

-'^ HI' —yfiiy^^ ap ^ -T fj,a y fJv 



(3.2b) 



T-imattcr 



diag(e,p_L,P^ P0,T^Prj) 



(3.4) 



Because of the conventions adapted in Eq. (3.1 ), the total 



energy and pressures are obtained by adding a factor of 
±i?2/2 to the matter terms. 



For the matter field, 



we first note that 13^ V'(^) 



S{x)V^4>{x)/^/t, where = d^, + U^ jx) + ieA^{x) 



is the covariant derivative defined in Eq. (2.23). So us 



ing the use the notation in App. \K\ th e field expansion 
(2.30), and the expectation value (2.46) of the mode op- 



erators, the diagonal components of the matter energy- 
momentum tensor (3.2a I are given by (no sum over /z) 



Tnmat 
MM 4 



k,k' A, A 

-1 



k,A 



(3.5) 



where we used the fact that Jfi{x) anficommutes with 
Il^{x), and the relation 'y'^Ilj^{x)j'^ — — n^(a;). Here we 
have defined the covariant derivatives 



V^, = d^,+Il^,{x) + ieA^{x) 
W^^d^,-n^{x) + ieAf,{x) 



(3.6) 



A. Energy Density 

For matter energy density term, Vq = dr and 7,- 
7"^ = 7", so using (|3.5|, we find 



sir) 



k,A 

^ETr[pL'^(r)i/k(r)] (3.7) 



2r 



fc.A 



So from ( |3.3| , the total energy density is given by 

^ = -E v,J-.)•pK/.K,-) + C• 



(3.8) 
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As it stands, the integral for the energy density in 
Eq. (3.8 1 diverges. We find the form of theses diver- 



gences by substituting the adiabatic expansion given in 
Eq. ( C20 1 and introducing a cutoff A,, in the tt,, inte- 
gral. This gives a adiabatic approximation to the energy 
density of 



E 



TTjj, the only terms that survive are 

^2 A 



kl 



2 

2 ,,,5 



P 



1 TT, 



1 



A 



12 V27rr 



A 

E( 



The first term renormalizes the field 

e^E^ El 



[l + e2(<Se2); 



2el 



(3.9) 



(3.10) 



The second term contributes to the cosmological con- 
stant, as we will see later. The third term is related to 
the zero-point energy of pairs of fermions. We regularize 
the energy density by computing the difference between 
£ and £a, £^^^ = £ ^ £a: which is now finite. 



B. Transverse pressure 

For the matter transverse pressure term, we have 7p 
-^P — —7^, and V p = dp, so from (3.51, we find: 



p^(^) = -^EM<^k'^^^)7Va;0L'^ (x)] . (3.11) 

So here we will need to find 

+00 

E xl^,m,hiP) 9p Xk^^m,-hiP) 

+ 00 

= h ^ Jrnikj_p) dp Jm+i{k±p) = hk±_ , 



m— — 00 



m— — 00 



where we have used the relation 

2 J,'„(z) = J„i-l(z) - Jm+l{z) 



Then (3.111 becomes 



k 



(3.12) 



and from (3.3 1, the total transverse pressure is given by 

E^ 



E 



hk 



(3.13) 



In order to study the divergences in the transverse pres- 
sure, we substitute the adiabatic expansion (C21b) into 
([3l3|. This gives 



Using Eqs. (2.651, due to the symmetric integration over Va-^ = 



E' 



A r 

E 



2uj 



kl 



1 TT^ 

"8 Cj5 



hk± M TTri 



1 TT,, TTjj 



(3.14) 



5 TT^TT^ 
8 



The second term in the above sum over p is odd in h and 
therefore vanishes. From Eqs. (2.65 1, the only terms that 
survive the tt,, integration are 



E' 



E 



e'E^ 



kl 


kl 






2 




1 1 


5 



1 / 3 71-2 5 TT,^ 



8 uj^ 8 
E^ 



(3.15) 



A 

E 



2lu 



From (3.10 1, the first term renormalizes the electric field. 
The term proportional to 1/r^ vanishes. Again, the 
transverse pressure is regularized by subtracting the adi- 
abatic expression from the divergent one, 7?™*^ — — 



C. Shear Pressure 

For the shear pressure term, we have 79 = — /9^7^ = 
— P7^ and Hg = 7^7^. The covariant derivatives (3.6 1 are 
given by 

7VV, = 7V5« + 7V/2, 
7VV, = 7V5«^7V/2, 

where 



cr^ 
Ut 



Uy 

(Jy 



So from (3.51, we find 
Pe{T) 



—I 
Apr 



;^ Aj^^L'^^x) [(7V9. + 7V/2) 0L^)(X)] 

[ ( 7°7^ do - 7°7 V2 ) ^l'^ {x) ] ' 01,^' [x) } . (3. 16) 
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So here we need to compute 

+00 



So terms proportional to 7^ vanish. Then from (3.5 1, we 
are left with 



m— — 00 



4-00 



{2m + l) Jra{k±p) Jm+l{kj_p) 



hk±p 



m — — 00 



where we have used the relation 

2mJra{z) = z[Jm+l{z) + Jm-l{z)\ 



Then (|3.16|) becomes 
1 hk 



P^(^) = 47^ E ^ {-^L'^ [7%' ( + eA{T) ) (x)] 

k,A 

+ [ ( + eAir) ) ^^'^ (x) ] ^ ^L"' } 
= 7^E[^"-^^(^)]^i^^(^) (3.20) 

= -E'^')^i;'='i''»('^')'^) • 



Adding the field pressure, we find for the total longitu- 
dinal pressure 



£;2 



(3.21) 



adding this to the shear pressure of the field, we find 

^2 



(3.18) 



Note that Ve = 'P±. The shear pressure is renormalized 
exactly like the transverse pressure. 



D. Longitudinal pressure 

For the longitudinal pressure {p — rj), A^{x) — —A(t 



and 7^ 
are given by 



-T7 . The covariant derivatives (3.6) 



« 7°7' = [ ^a,^ + eA{T) ] + 7V2 , 
z 7%' = 7"73 [ + eA(r) ] - 7V2 , 



where 



„0„,1„2„,3 



1 



(T, 



7 =i7 777 =1^^1, 77 _ ^ 



Here 7^ flips the upper and lower components of the 
spinor, which leads to the equation 



+00 



m— — 00 



+00 

5 E [ j'^^f'^P) - Jrn+lik±p) ] = . (3.19) 



From Eq. (|C21a|), the adiabatic expansion of the longi- 

(3.22) 



tudinal pressure is given by 
_V[^-(fci + Af2)(i^-^^) + 



From Eqs. (2.651, the only terms that survive the tt^ 
integration are 



V 



2 ^ ^7:1 5 



E 



V 



e^E^kl + M^)^ 

UJ' 



p 

„2 I fl ,r2 ■ 



2r2 





5 




V w5 


" 4 





(3.23) 



Again, the first term renormalizes the electric field and 
the second term renormalizes the cosmological constant. 
The finite part of the longitudinal pressure is given by 
"^r*" = "P,, - Va;ri, as before. 



E. Conservation equations 

The covariant derivative of the energy-momentum ten- 
sor in boost-invariant coordinates is conserved 

T'^'.f, = a^T^" + Ti^^T"" + Tf^^T^"" = . (3.24) 

The only nonzero Christoffel symbols are given in 
Eq. (2.251. There are only two conservation equations 



that result from Eq. (3.24). For — t, (3.24) reduces to 



or dr (tS) +Vr,^0. (3.25) 
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Using the equation of motion (2.39 



tion (2.63 1, one can show that Eq. (3.25 1 is automatically 



and Maxwell's equa- 



satisfied. 

The conservation equation for v — p amounts to a 
relation between the transverse and shear pressures. We 
find that V \_ = T'e, which is satisfied by our expression 
in Eqs. (3.131 and (3.181. For the one-dimensional boost 



invariant expansion we had instead for the equation for 
the energy density, 

{tE) + 7^ = . (3.26) 



IV. QUASIPARTICLE PHASE SPACE 
DISTRIBUTION FUNCTIONS 

For the problem at hand, particle production from clas- 
sical electric fields, it is possible to introduce an inter- 
polating number density via a Bogoliubov transforma- 
tion that is an adiabatic invariant. This was done previ- 
ously [HJ [551 1231 [21] )• However, when we consider the 
fully interacting case with quantum gauge fields then one 
needs to resort to an "effective" quasiparticle distribution 
function that allows one to reproduce the expectation 
value of the current and the energy- momentum tensor. 
That is, we want to determine an effective distribution 
function /(x, fc) in analogy with relativistic kinetic the- 
ory (see for example Refs. [571 [IH]) such that 



( J^) = e / mk''f{x,k) 



and 



where 



d^k 



Dfc = 2r6'(fc°)(5(fc2 - M^) 



(27r)3 



(4.1) 



(4.2) 



(4.3) 



Here, r is a degeneracy factor which counts the number of 
species. For our case of quark and anti-quark pairs with 
spin one-half, we have r — A. Hence, the renormalized 
comoving energy density, £ — Too, is given by 



+00 



kj_ dk± 
2n 



+ 00 



dfc. 



- ^k^,k^{T) f{T,kj_,kr,) 



(4.4) 

and will be identified with the renormalized field theory 
result: 



2 



, (4.5) 



where we have subtracted the divergences coming from 
the cosmological term and the charge renormalization. 
Note that when the single-particle distribution becomes 
independent of proper time, r, one can easily derive the 
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FIG. 1; (Color online) Convergence of the proper-time evolu- 
tion of the electromagnetic fields, current, energy and lon- 
gitudinal pressure with respect to the cutoff A, in (3-1-1)- 
dimensional QED. Here we choose m = 1, A{tq) = and 
E{to) = 4. 



conservation of energy equation in terms of the energy 
density and longitudinal pressure: Consider the identity 



4 /■+°°fc_Ldfc_L f+°°dk 



2tt 



r, df{T, k^,k„) 
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FIG. 2: (Color online) Early proper-time evolution of the 
electromagnetic fields and current with respect to the cutoff A, 
in (3-|-l)-dimensional QED. We note that for r < 5 the results 
are dependent on the cutoff, whereas at later proper times the 
results become insensitive to the cutoff. This is an artifact of 
our choice of initial conditions, which are not consistent at 
early times with the adiabatic-expansion-based substraction 
scheme. 



then integrate by parts to obtain 

OT T 

where the longitudinal pressure is introduced as 



Vn 



4 



27r 



27r LOk^^k^ir) 



(4.7) 

The quasiparticle phase-space distribution of pairs of 
particles and antiparticles with a specific spin in light- 
cone variables is introduced as 



(Pxj_dri d'^kj_dk^ (27r)3 
such that the pair density is obtained as 
d^iV /•+°°fc_Ld/cj f+°°dk 



(4.8) 



d'^x^di] 



2n 



2tt 



^f{T,k^,k,). (4.9) 



Here, we have d^x± = pdpdd and d^k± = k±dk±d(j). 
For completeness, we note that in (1+1) dimensions the 
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FIG. 3: (Golor online) Proper-time evolution of the transverse 
(fex-projected) distribution of the current for r values of 2, 3 
and 9. The oscillations of the transverse distributions of the 
current present at early proper times are consequences of the 
inconsistency between the adiabatic-expansion regularization 
scheme and our choice of initial conditions which are designed 
to keep the storage and time requirements of the simulation to 
a minimum. At later proper times, these oscillations dampen 
out and eventually disappear. Therefore, the early proper- 
time evolution should be regarded as "unphysical" and will 
be disregarded. 



pair density reads 

diV 
dr] 



-t-oo 



dk 



27T 



^ fir, fc„) 



(4.10) 



When the pair distribution becomes independent of the 
proper time r then we are in the "out-regime" and can 
stop our calculation as far as determining the particle 
spectra. We need to relate this quantity to the center-of- 
mass distribution of electrons and positrons produced by 
the strong electric field. We introduce the free-particle 
rapidity, y — ^ln[(i? + kz)/{E — kz)], and "transverse" 
mass, M±_ = \/k\ + m?, by relating them to the carte- 
sian coordinate four-momentum in the center-of-mass 
system, A;° = (-E-jk), by the relation 

k°- = ( M_L cosh y , k_L , Af_L sinh y ) 
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FIG. 4: (Color online) Proper-time evolution of the elec- 
tromagnetic fields and current for the case of (1-1-1)- and 
(3-|-l)-dimensional QED, respectively. Here we choose m = 1, 
A{to) = and E{to) = 4. 



The boost that takes one from the center-of-mass coordi- 
nates to the comoving frame, where the energy momen- 
tum tensor is diagonal, is given by tanhry = v = z/t, 
so that one can define the "fluid" four-velocity in the 
center-of-mass frame as 



= (cosh T], 0, 0, sinh ry) 



(4.11) 



It is important to relate the momenta canonical to rj and 
T to the center-of-mass variables. In the out regime we 
can identify these canonical momenta from the free parti- 
cle (l-l-l)-dimensional Lagrangian in covariant form. We 
show now that 



1 



■In 



t + z 



(4.12) 



have as their canonical momenta 
/ct- — Et/r ~ kz z/t , 



k,^ = -Ez + tkz. (4.13) 



Consider the metric ds^ = dr^ — r^drf' and the free par- 
ticle Lagrangian in (1+1) dimensions 



(4.14) 
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FIG. 5: (Color online) Proper-time evolution of the matter 
components of the renormalized energy-momentum tensor for 
the case of (1+1)- and (3+l)-dimensional QED, respectively. 



Then we obtain 
dr 

kr^M M 
as 





\ dt 




dz 






[dz, 





(4.15) 



{Et — kz z) — k'^Ua = M± cosh(ry - y) , 



and 

krt 



ds 



dt 
ds 



drj 



dz ) ^ ds 



(4.16) 



'E z -\- kzt ~ — r Af_L sinh(77 — y) 



It follows that kr ~ A;^ = M±_ cosh(77 — y) has the 
meaning of the energy of the particle in the comoving 
frame. 

The interpolating phase-space density / of particles de- 
pends on (t, k±,kri) and is found to be ?7-independent. In 
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FIG. 6: (Color online) Proper-time evolution of the den- 
sity of pairs, dN/dy, and proper time evolution of the ratio 
Te{T)/[dN/d'y] for the case of (1-fl)- and (3-|-l)-dimensional 
QED, respectively. The circles in the upper panel denote step 
increases in the particle density corresponding to the second 
and third gradients of the current. 



order to obtain the center-of-mass particle rapidity and 
transverse momentum distribution, we change variables 
from fc^) to {z,y) at a fixed proper time r, i.e. 



d^iV 



d'^x±(Pk±dz dy d'^x±d^k±dr] dk,i 
we have 



d{z,y) 



So, from (4 
d^N 



(4.17) 



d'^k±dy 



d^x± dz 



djy, kr,) 
d{z,y) 



f{T, kj_,kr,) 
(27r)3 



(4.18) 

where the Jacobian is evaluated at a fixed proper time r. 



(4.19) 



d{r],k^) 




dkjjdy dk,f/dz 




dk,f drj 


d{z,y) 




drj/dy drj/dz 


T 


dy dz 



M±^ cosh(?7 — y) 
cosh rj 

However, since at fixed r, we have 



dk,^ 




dkn 


dy 


T 


drj 



we obtain 



d{r],kn) 


dkjj 


diz,y) 


^ dz 



(4.20) 



(4.21) 



Calling the integration over the transverse dimensions the 
effective transverse size of the colliding ions A± = TrR^g 



1 f 

J dkrjf{T,kj_,k^) 



d^N 
d^fcj^ drj 



we then find from (4.181 that: 

d^N _ A_ 

d^kZfhj ~ J2n) 

(4.22) 

The quantity in Eq. (4.221 is independent of y which is 
a consequence of the assumed boost invariance. There- 
fore, using the property of the Jacobean, we have proven 
that the distribution of particles in particle rapidity is 
the same as the distribution of particles in fluid rapid- 
ity, verifying that in the boost-invariant regime Landau's 
intuition was correct |29j . 

We now want to motivate the Cooper-Frye formula 
used to calculate particle spectrum in hydrodynamical 
models of particle production [27] . We have that a con- 
stant T surface, which is the freeze-out surface of Landau, 
is parametrized as 



dE° = (dz,0,0,dO (4.23) 
= Aj_ T drj ( cosh 77, 0, 0, sinh rj ) . 

Therefore, we find 

k'^dY.a^ Aj_Mj_T cosh{T]-y)dr] (4.24) 
= A± I dk„ I . 



Thus, we can rewrite our expression for the field theory 
particle spectra as 

d^iV 
d^fcj^ dy 



(27r)3 



dkrf f{T, kj_, kr,) 
rdS,/(T,fci,^) , 



(4.25) 



where in the integration we keep y and r fixed. Thus, 
with the replacement of the thermal single-particle dis- 
tribution by the quasiparticle distribution function, we 
get via the coordinate transformation to the center-of- 
mass frame the Cooper-Frye formula. For completeness, 
we note that in (1+1) dimensions the particle spectra, 
dN/dy, are given by the integral in Eq. ( 4.10[ ). 

The boost invariant assumption leads to an energy mo- 
mentum tensor which is diagonal in the (r, p, 0, 77) coor- 
dinate system which is thus a comoving one. In that 
system one has for the matter energy-momentum tensor 



(4.26) 



Thus we find in this approximation that there are two 
separate pressures, one in the longitudinal direction and 
one in the transverse direction which is quite different 
from the thermal equilibrium case. However only the 
longitudinal pressure enters into the energy conservation 
equation: 



d(rf) 



E J„ 



(4.27) 



It is useful to rewrite the conservation of energy in the 
out regime as : 

d£ £ + V„ 
dr T 



. 



(4.28) 
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So to the extent that the ultra-relativistic one 
dimensional equation of state £ 
has the simple result 



Vn is true, then one 



E 



(4.29) 



It turns out, as our numerical results show below, that 
although pri/e ~ 1 for part of the period of the oscillation, 
during the minima — s- and this seems to be governing 
the fallofF which is more like £ oc 1 /r. 



V. NUMERICAL RESULTS 

Let us review the equations we intend to solve numeri- 
cally. The first of these is the polarization equation ( 2.39 ) 



9.Pi')(r) = 2k,(r)xPi')(r), 



and the second is the backreaction equation (2.71 ) 

A r 



drE{T) 



l + e2((5e2) 
eE{T) 



E 



(5.1) 



(5.2) 



'^ + ^(fci + M^)(-L-'4) 
uj T ^ ^ ' \ 4^5 4lj7 / 



Here we have subtracted from the integral the adiabatic 
expansion of P['^]}^ (tt^ , r) . 

In order to solve the coupled Dirac and backreaction 
equations, we construct a grid in k± and k^j space as 
follows: The fc^-momentum variable is discretized on a 
nonuniform piece-wise momentum grid with a cutoff at 
kjj — Arf] we find that a value of « 500 is necessary to 
obtain numerical results insensitive with respect to the 
cutoff. A similar nonuniform grid is used to discretize the 
k± variable, fcj_ g [0,A]. A fourth-order Runge-Kutta 
method is employed to solve the coupled Dirac equation 
and backreaction problem. 

For the purpose of calculating the subtracted values 
of the current J(t) and the components of the matter 
energy-momentum tensor, we compute the momentum 
integrals symmetrically with respect to the variable tt^ 
rather than /c^. The corresponding momentum cutoff in 
7r,,-space is chosen to be 20% greater than TmaxA^ to 
allow for possible very large values of A{t), the latter 
being unknown at the beginning of the calculation. 

The momentum integrals with respect to tt^ and k± are 
performed using a Chebyshev integration method with 
spectral convergence [31 . Using this procedure, we found 
that a grid of approximately 8000 points in the fc^ (or 
TTri) variable and 128 grid points in the k± variable is 
necessary to obtain a converged numerical result. As 
such, the calculations for the backreaction problem in 
(3-fl) -dimensional QED require at least 100 times larger 
storage and computational time then the corresponding 
(l-l-l)-dimensional QED problem. 

For illustrative purposes, we took: m = 1, e = 1, 
To = 1/m = 1, A{to) = 0, and E{to) = 4. These 



strong-field initial conditions have been shown to pro- 
duce sufficient fermion pairs at r = tq for plasma oscil- 
lations to take place. Just like in the (l-l-l)-dimensional 
case, the conservation of the energy-momentum tensor, 
see Eq. (3.251, serves as a numerical test: for the results 
of simulations reported here, the renormalized energy- 
momentum tensor is conserved within machine precision. 

In order to keep the size of the simulation to a min- 
imum, we chose the initial conditions corresponding to 
the one-field scenario introduced first in Ref. |24| and 
summarized in Sec. |IIC"| In Figs. [T] we illustrate the con- 
vergence of our results with respect to the choice of the 
cutoff, A. For completeness, we depict the proper-time 
evolution of the fields, A{t) and E{t), current, J(t), en- 
ergy, e(r) and transverse pressure, Pri^T), for cutoff values 
between 2 and 6. We conclude that for A = 5 the results 
are insensitive to the cutoff A, within numerical accuracy. 

It is important to note that our choice of initial con- 
ditions results in a time evolution that is not consis- 
tent with the adiabatic expansion for early values of the 
proper time. In Fig. [2] we depict the proper-time evolu- 
tion of the fields, A{t) and E{t), and current, J(t), at 
early times for several values of the cutoff A. We note 
that while the proper-time dynamics converges for r > 5, 
for earlier times the proper-time evolution depends on the 
choice of the cutoff A. However, for larger proper-time 
values, the nonadiabatic components of the current dissi- 
pate and the adiabatic-expansion-based subtraction be- 
comes exact. This behavior is illustrated numerically in 
Fig. [3j where we depict the fc_L-projected distribution of 
the current for r values of 2, 3 and 9. We notice that the 
nonadiabatic oscillations of the current present at early 
proper times dampen out and disappear at later proper 
times. Therefore, the early proper-time evolution will be 
disregarded as "unphysical." This is a small price to pay 
in order to keep the storage and time requirements of our 
simulation to a minimum. 

In Fig. g] we compare the proper-time evolution of 
the electromagnetic field, A{t), electric field, E{t), and 
current, J(t), for (1+1)- and (3-|-l)-dimensional QED, 
respectively. Similarly, in Fig. [5] we depict the proper- 
time evolution of the matter components of the energy- 
momentum tensor. We note that in (1-1-1) dimensions the 
fields evolve much faster than in (3-1-1) dimensions. Also, 
in (3+l)-dimensional QED the energy and longitudinal- 
pressure densities are very very close in magnitude, which 
in turn results in a small transverse pressure, p±(t). 
Qualitatively, we also note that the modulation observed 
in the proper-time evolution of the current and longitu- 
dinal pressure in (l-l-l)-dimensional QED are not present 
any longer in (3-1-1) dimensions. Also by inspecting the 
two upper panels in Fig. [5] we notice that in 3 -I- 1 di- 
mensions Pri ~ e. However, the ratio Prj/e becomes close 
to zero near the minimum of the oscillation and this in 
turn leads to re to be almost constant instead of going 
as 1 /r using the arguments coming from the energy con- 
servation equation Eq. (4.281. 

Finally, the proper-time evolution of the density of 
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FIG. 7: (Color online) Proper-time evolution of the longitudinal momentum-dependent pair-density distribution, n^^ , and the 
proper-time evolution of the transverse momentum-dependent pair-density distribution, n^j^ . (See also Ref. I30p 



pairs, diV/dy, are depicted in Fig. [6] Particles are be- 
ing created corresponding to the current gradients, with 
the major contribution corresponding to the initial cur- 
rent gradient, and subsequent smaller step increases be- 
fore the particle density saturates. At late values of the 
proper time, the ratios re(r) /[dA^/dy] are seen to ap- 
proach a constant consistent with the hydrodynamical 
picture, which relates the energy in a bin of rapidity di- 
vided by the energy of a single particle with that rapidity 



with the number of particles in a bin of rapidity as ex- 
plained in Ref. 23. In the real problem we expect that 
interactions between the fermions will eliminate the os- 
cillations observed here. 

The proper-time evolution of the momentum- 
dependent longitudinal pair-density distribution, Ut^^ , 
defined as 

"-^^^dl^^i ^/(r,fc^,M, (5.3) 



16 




FIG. 8: (Color online) Comparison of the transverse distri- 
bution of particles according to our numerical simulation and 
the transverse distribution of particles for the constant elec- 
tromagnetic field case with eE = 4 as described by Eq. ( |1.2[ ). 



dimensions suggesting that the extra degrees of freedom 
perform some smoothing. We now have the first nu- 
merical results for the transverse momentum distribution 
function of fermion pairs which we can compare with the 
exact results for the constant field problem. We find 
that unlike the constant field case, the distribution is bi- 
modal. At modest fc^ < 5m^ the transverse distribution 
is similar to the constant field case with a reduced (75%) 
effective eE for cEq = 4. For larger transverse momen- 
tum k'j_ > lOm^ the transverse distribution function has 
a tail described by an effective eE which is of the order 
of 50. This is a totally new feature that is as of yet not 
understood simply. In a related paper |32j we will also 
consider a transport approach to the (3-|-l)-dimensional 
problem and show that such a semiclassical picture works 
better in (3-1-1) than in (1+1) dimensions. 
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and the transverse pair-density distribution, n^^ , defined 
as 



A, 



d^A^ 
d^fcj^ dy 



+ C30 



dfc,j 



f{T,k^,k^), (5.4) 



in (3+l)-dimensional QED are shown in Fig. [Tj (See 
also Ref. 30 ) We note that the centroid of the particle- 
density distribution, n^r^ , oscillates between positive and 
negative values of tt,,, similar to the (l-l-l)-dimensional 
QED case. 

In Fig.|8]we compare the transverse momentum distri- 
bution given by the constant field exact solution Eq. (1.2) 
with the results of our numerical solution for eE = 4. 
Part of the results are expected, that is at small trans- 
verse momenta the distribution of particles is similar to 
the static case but with a smaller effective field since the 
field is decreasing during the first phase of particle pro- 
duction. What is unexpected is that, in the problem with 
backreaction, there is a new tail in the transverse momen- 
tum distribution which falls exponentially with an effec- 
tive \eE\ = 50. This is a totally surprising result whose 
origin we do not yet have a simple explanation for. 



VI. CONCLUSIONS 

We have for the first time calculated the transverse 
distribution of jets produced by an initial strong elec- 
tric field including the effects of backreaction. We have 
compared the results of our (3-1-1) -dimensional calcula- 
tions (for "hydrodynamic" quantities as well as for the 
proper time evolution of the electric field and current) 
with their (l-l-l)-dimensional counterparts. We find that 
the electric field degrades much quicker in (3-1-1) dimen- 
sions than in (1+1) dimensions. Also secondary oscil- 
lations in the current and in the longitudinal pressure, 
present in (1+1) dimensions seem to be absent in (3+1) 
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the United States Department of Energy. The authors 
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APPENDIX A: NOTATION 

In this appendix, we list our notation and conventions 
used throughout this paper. We use a boldface k to desig- 
nate the complete set of mode variables, k and k,n subsets 
of the full set, and k^ the cylindrical coordinate set. In 
addition, the set p substitutes the tt^ kinetic momentum 
for kri- These sets are given by 

k = ( fcj^, /i, /c^, m ) , (Ala) 

k={k^,h,kn) , (Alb) 

p={k±,h,7Tr^), (Ale) 

km = {k±,m), (Aid) 

k(p = ( kxj ky^k^^ — ( kj_^ (j}^ k^i ) . (Ale) 

Sums over these quantities indicate the following inte- 
grals and sums 



E 

k 

E 

k 

E 



E 



+ 00 



k_\_ dk±^ 
2n 



+00 



+00 



27r ^ ^ 

m= — oo h=±l 



27r 



+ 00 



dfc, 



} E 



2^ 

+°°dfc, 
-oo 27r 

°° k\ dk 



h=±l 

+°°d7r^ 



+00 



27r 



dky 



2-K 



n 



d^ 

-oo 27r 

+°°dk^ 
27r 



(A2a) 
(A2b) 
(A2c) 
(A2d) 
(A2e) 
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We also use the following notation for 5-functions 



^k,k' 



^6{k,-k'^). (A3) 



In a similar way, we put x = ( p,9,ri). The sum over x 
means 



E 



pdp r^^'de r+°°dr) 

^ Jo 2n 2n 



(A4) 



and (Sx.x' means 



6^,^,^^-^^^S{0^0')S{v-v')- (A5) 



APPENDIX B: TRANSVERSE HELICITY 
EIGENVECTORS 

In this section we derive transverse helicity eigenvec- 
tors and show how they can be used to expand solu- 
tions of the Dirac equation in boost-invariant coordi- 
nates. The Hermitian two-component transverse helicity 
operator H± is defined in momentum space by 

where — k± cos (f) and ky — k± sin (j). We write the 
eigenvalue equation for this operator as: 

H±X<t,,h ^ hx<t,M , (B2) 
with eigenvalues h = ±1 and orthogonal eigenvectors: 

X<p,h = ^ ■ (B3) 

We note that az X4>,h — X<p-h, and that: 

{(Txkx + cryky)x4,.h^-ihk±x<p-h- (B4) 

For the coordinate system x'^ = {T,x,y,ri), Dirac 's 
equation is: 

{ rir, v) [td^-gA^iT)]-M}4,iT,x,y,v)^0, (B5) 
where now: 

I'^iv) — cosh rj^^ — sinh rj 7'^ , 
7'' ('''i ^) = ( ^ sinh ^77'' + cosh rj'y^^ /t . 

with 7^ = 7^ and 7^ = 7^ . We next boost to a coordinate 
system where 77 = by setting: 



'4'{T,X,y,T]) = S{r]) (p'{T,X,y,T])/y 

then (|B5|) becomes 



(B6) 



+ 7' [ i + A{t) ]/t-M} <^'(r, a;, y, 77) = , (B7) 



which is what we want to solve. So let us first introduce 
the Fourier transform: 

(/.'(r, X, y, 77) = e^(''^-+'^yy+'^^^) , (B8) 



where k = (kri,kx,ky). Then (B7) becomes: 
[ * 7° 9, - 7I fc, - 72 ky ~ 73 (r) - M ] (^'fc (r) = , (B9) 

where 7r^(T) = [fc^ — eA(r) ]/t is the kinetic momentum. 
Using Eq. ( B4 1 , we see that ( B9 1 is separable if we put 



<j^',{T) ^ ('^j+y-'^\^\^^-+') , (BIO) 

where 0(±);fe(T) now satisfy the two-component equation 



Tr^{T) + ihk^_ -M J l^(/)(_);fe(r) 



which agrees with Eq. (2.351. Near r — tq, there are 



positive and negative energy solutions to these equations 
which we label by A = ±1. So the fermi field 0(r, x, y, rj) 
can be expanded as 



(A) 



h=±l X=±l 



'f'l+yA^) X4>, + h \ j( k^,^+k^ x+ky y ) 



(Bll) 



where A^^} h ^'^^ creation and annihilation operators 
for the state described by (fc^, h, A). Now let us introduce 
cylindrical coordinates, x = p cos 9 and y = p sin 9, so 
that 



kx X + ky y = k±p cos(0 - 



(B12) 



Now the generating function for Bessel functions is given 

by 



+00 



exp[z(i-lA)/2]- Y t"'J,n{z) 



(B13) 



m— — 00 



If we put t — ie^^ and z — k± p, this becomes 

4- 00 

exp[z/fcj^pcos(6'-0)] = Y ^'"e™^^-'^) J™(fc_Lp) 



Using these results in Eq. (Bill, we find in cylindrical 
coordinates the expansion. 



+00 



i'ir,pAv) = Y E E E4ie-- 



kj, m= — 00 /i=±l A=±l 



0|?)^,(r)x0,+a^™^,„,(,.,)^^^(^^^^^ (B14) 

'/'(_);fe(T)X0,-V 
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Now let us define the Fourier transform pair: 

e2TT 



(A) 

,k± ,m,h 



A 



(A) 







2^ kr,,k±,<ji,h 



kj^ ,k± ,m,h 



-inKp 



(B15a) 
(B15b) 



m— — OO 



and fourth components, Eq. (|B14|) becomes 



k A=±l 



X e' 



where now 



(B16) 



(B17) 



Finally, we boost to a coordinate system where = by 
multiplying 0'(t, p, 6, rf) by Sp^{9), which can be written 
as 



S~^{9) = cos(6l/2) - sin(6i/2) 
















g-ie/2 














g+ie/2 





V 











(B18) 



so that ^{x) — S{9,rf) 4>{x) and from Eq. (B16), we find 

(B19) 



,)^s-\e)4>'{x) = Y, E ^ 

k A=±l 



where 



and where 



„i(m+l/2)e 



~h Jjn+l{kj_p) 



(B20) 



(B21) 



in agreement with the field expansion given in 
Eqs. (|2.30|, (|2.32|, and (|2.33|) in Sec. Ill Bl We have shown 



here that the separation of variables method for the Dirac 
equation we used in Sec. |IIB| can easily be understood as 
an expansion of transverse helicity eigenvectors in boost- 
invariant coordinates. 



The transverse helicity eigenvectors given in Eq. (B21 ) 
satisfy the eigenvalue equation, 



de 



-'hk^Xk^,h{p,e), (B22) 



are normalized. 



00 1^27: 

pdp / d9Xk„,,hiP^0)xk' ,h'iP^^) 
a Jo " 

6{k^ - k'J 



(B23) 



So using (B15al, and putting m ^ m + 1 in the second and complete 



2tt 



-\-oo 



E E x,,^Jp,o)xLAp',o') 

,n S{P-P') 



m= — oo h=±l 



6(9^9') 



'PP 



(B24) 



APPENDIX C: ADIABATIC EXPANSION OF 
SOLUTIONS OF THE DIRAC EQUATION 

In this section, we find an adiabatic expansion of the 
positive energy solutions of the Dirac equation for a 
slowly varying field A{t). It is simplest to obtain an 
adiabatic expansion of the polarization vector P^^^(r), 
which we introduced in Sec. |IIB[ The equation of mo- 
tion of the polarization vector was given in Eq. (2.391 
as 



Pi'Hr) = 2k,(r)xP^^^(r), 



(CI) 



where k/j(r) is given by 

kfe(r) = 7r^(r)ei-|-/ifc_Le2 + Afe3. (C2) 



The initial condition at t = tq is given in Eq. (2.531 as 



(C3) 



For slowly varying values of 7r^(T), 'P^^\t) simply pre- 
cesses about the slowly varying value of kfc(T). In order 
to count derivatives with respect to t, let us put 



We next expand P^^-'(t) in powers of e by writing 

P = p(0) + eP(i) + e2p(2) + ... _ 
P = eP(°) + e2p(i) + e3p(2) + ... , 



(C4) 



(C5) 



Here and in the following, we omit momentum and time 
dependencies and the (A) label. The superscript now 
counts powers of e and the dot refers to derivatives with 
respect to t. So substitution of (C5) into Eq. (CI I be- 
comes 



ep(0)+,2p(l)^^3p(2)_^... 

= 2kx [p(0) + eP(i) + e2p(2) 



(C6) 
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Equating equal powers of e gives 



2k X P(°) 


= 0, 




(C7a) 


2k X P(i) 


= P(0) , 




(C7b) 


2k X P(2) 


- -6(1) 

-■- J 




(C7c) 


2k X P(3) 


- p(2) 


etc • • • 


(C7d) 



Let us introduce transverse and longitudinal components 



L ' 



of the polarization vector by writing P*-"' — j. 
where k • P^"^ = and k x P^''^ = 0. So Eqs. 



determine only the transverse components of the polar- 
ization vector. The longitudinal portion is then fixed by 
the normalization requirement, as we will see below. 
From Eq. (C7al, P'"' is entirely longitudinal and has 



the normalized solution 

P(°) = k/w , L0= |k| = 

So P('') = Pq. So we have 



kl 



p(o) ^ k _ k^ 
uj up- 



k k (k • k) k X (k X k) 



Then Eq. ( |C7b| becomes 

^^^\ k X (k X k) 
2k X P(i) = , 

SO the transverse component of P'^-* is given by 
t,(i) _ k X k 



2^3 



We will choose the longitudinal component P]^ 
that 



(1) 



(C8) 
, (C9) 
(CIO) 

(Cll) 
= 0, so 



p(i) 



k X k 

2cj3 



TT, 
2UJ 



"^{Me^-hk^e^) . (C12) 



Then to first order, P = p(°) + eP^^), the polarization 
vector is normalized to this order, since 



p2 = p2 + 2eP(°) -P^i' = 1 



(C13) 



From (C12l, we find 



p(i^ k X k 3 (k X k) cO 

_ k X k 3 (k X k) (k • k) 
3(k-k)k-u;2k 



(C14) 



k X 



2^5 



Then Eq. ( C7c ) becomes 
2k X P(2) ^ k X 



3 (k • k) k - w2 k 

2^ 



(C15) 



so adding a longitudinal part to P(^\ we find 

p(2)^3(k-k)k-C.2k 



'AAok. 



(C16) 



where Af2 is to be fixed by the normalization requirement. 
From the expansion ( C6 ) , we find to second order 



P2 = l + e2 [p(i)2^2PW -P^^)] 



1 . (C17) 



So we want to choose 7V2 such that p(i) 2 + 2 P(°) • P^^) ^ 
0. This gives the equation 



p(o) . p(2) ^ 3(k-k)2-a>2k.k 



-M2UJ 



2 8cj6 ' 

from which we find 

3(k-k)2-t^2j^.k |k|2cj2 _ (k.k)2 



(C18) 



Af2 = -- 



1 < 



1 TT^ TTrf 



c ^2 -2 

5 7r„7r„ 



(C19) 



8 ^5 4 ^5 8 ^7 ■ 
So to second adiabatic order, the polarization vector is 

given by 

p = p(0)+,p(l)+g2p(2)^... (C20) 



k 



k X k 



3 (k • k) k - k 



2u;3 

In component form, we find 



4^5 



'M2k 



w ^ ^ V 4 

hk± 



5 TT^ TT^ 



8 OJ^ 



P 



2a;3 



M , , TT, 

UJ 2 Cj3 

■2 



1 tt: 

8 cj5 ' 4 uj^ 



1 TT,, 7f,, 5 71,2 



V 8 cj5 



1 TT^ TT^ 

4 Cj5 



8 UJ^ 



5 TT^TT^ 
8 



(C21a) 
(C21b) 

(C21c) 



which completes the adiabatic analysis used in this paper. 



[1] T. S. Biro, H. B. Niclson, and J. Knoll, Nuc. Phys. B245, 

449 (1984). 

F. Sauter, Z. Physik 69, 742 (1931). 
W. Heisenberg and H. Euler, Zeit. fiir Physik 98, 714 

(1936). 

J. Schwingcr, Pliys. Rev. 82, 664 (1951). 
B. Audcrsson, G. Gustafson, G. Ingelman, and T. Sjos- 
trand, Phys. Kept. 97, 31 (1983). 

X.-N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 
(1991). 

L. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 
(1994). 

L. McLerran and R. Venugopalan, Phys. Rev. D 49, 3352 
(1994). 

D. Kharzeev, E. Levin, and K. Tuchin, Phys. Rev. C 75, 
044903 (2007). 

A. Bialas and W. Czyzn, Phys. Rev. D 30, 2371 (1984). 
A. Bialas, W. Czyzn, A. Dyrek, and W. Florkowski, Nuc. 
Phys. B296, 611 (1988). 

K. Kajantic and T. Matsui, Phys. Lett. 164B, 373 
(1985). 

G. Gatoff, A. K. Kerman, and T. Matsui, Phys. Rev. D 
36, 114 (1987). 

V. V. Skokov and P. Levai, Phys. Rev. D 78, 054004 

(2008) . 

G. C. Nayak and P. van Nieuwenhuizen, Phys. Rev. D 
71, 125001 (2005). 

G. C. Nayak, Phys. Rev. D 72, 125010 (2005). 
M. Claudson, A. Yildiz, and P. H. Cox, Phys. Rev. D 22, 
2022 (1980). 

A. Cashcr, H. Neuberger, and S. Nussinov, Phys. Rev. D 

20, 179 (1979). 

G. C. Nayak, Phys. Rev. D 72, 125010 (2005). 
F. Cooper and E. Mottola, Phys. Rev. D 40, 456 (1989). 
Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, and 

E. Mottola, Phys. Rev. Lett. 67, 2427 (1991). 
Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, and 

E. Mottola, Phys. Rev. D 45, 4659 (1992). 

F. Cooper, J. M. Eisenberg, Y. Kluger, E. Mottola, and 

B. Svetitsky, Phys. Rev. D 48, 190 (1993). 
B. Mihaila, J. F. Dawson, and F. Cooper, Phys. Rev. D 
78, 116017 (2008). 

F. Cooper, S. Habib, Y. Kluger, E. Mottola, J. P. Paz, 
and P. R. Anderson, Phys. Rev. D 50, 2848 (1994). 
B. Mihaila, J. F. Dawson, and F. Cooper, Phys. Rev. D 
74, 036006 (2006). 

F. Cooper, G. Frye, and E. Schonberg, Phys. Rev. D 11, 

192 (1975). 

E. A. Calzetta and B.-L. Hu, Nonequilibrium quantum 
field theory (Cambridge University Press, Cambridge, 
2008). 

L. D. Landau, Izv. Akad. Nauk SSSR 17, 51 (1953). 
Animations of the proper-time evolution of the 
longitudinal and transverse momentum-dependent 
particle-density distributions can be found at 
http://einstein.unh.edu/BogdanMihaila/qed-3-l-l/. 
B. Mihaila and L Mihaila, J. Phys. A: Math. Gen. 35, 
731 (2002). 

J. F. Dawson, B. Mihaila, and F. Cooper, in preparation 

(2009) . 



